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The static vacuum spherically symmetric solutions in massive gravity are obtained both analytically and numerically. 
The solutions depend on two parameters (integration constants): the mass M (or, equivalently, the Schwarzschild radius), 
and an additional parameter, the "scalar charge" S. At zero value of S and positive mass the standard Schwarzschild 
black hole solutions are recovered. Depending on the parameters of the model and the signs of M and S, the solutions 
may or may not have horizon. Those with the horizon describe modified black holes provided they are stable against 
small perturbations. In the analytically solvable example, the modified black hole solutions may have both attractive 
and repulsive (anti-gravitating) behavior at large distances. At intermediate distances the gravitational potential of a 
D modified black hole may mimics the presence of dark matter. Modified black hole solutions are also found numerically 
in more realistic massive gravity models which are attractors of the cosmological evolution. 
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1. Introduction 

Recent advances in observational cosmology [1, 2, 3] 
have revived interest in alternative theories of gravity in 
which the gravitational interaction is modified in the in- 
frared domain and which could — potentially — explain 
the accelerated expansion of the Universe without intro- 
ducing the dark energy and matter components. Theo- 
retical consistency and existing experimental data impose 
severe constraints on such models. Different approaches 
to the problem have been discussed in the literature (see, 
e.g., Refs. [4, 5, 6, 7, 8, 9, 10, 11], and Rcfs. [12, 13, 14] for 
reviews). One of them employs spontaneous breaking of 
Lorentz symmetry by space-time dependent condensates of 
scalar fields [15, IG] coupled to gravity in a covariant way 
via a derivative coupling. The resulting theory may have a 
non-pathological perturbative behavior about the broken 
vacuum [17, IG] and exhibit modifications of gravitational 
interactions at large scales (see Sect. 2 for more details). 
In particular, graviton may acquire a non-zero mass, which 
is the reason to call these models massive gravity models. 

Existing experimental data constrain the mass of the 
graviton and other parameters of massive gravity mod- 
els. Absence of Lorentz invariance makes the constraints 
weaker than one would expect in a Lorentz-invariant the- 
ory: the Newton's potential remains unmodified in the lin- 
ear approximation despite non-zero graviton mass [18], so 
that Solar system constraints are satisfied for rather large 
masses. The constraint on the mass of the graviton comes 
from the emission of gravitational waves by binary pulsars 



which is responsible for their spin-down [19]. Consistency 
of observations with GR requires the graviton mass to be 
smaller than the inverse period of orbital motion of the 
binary system, that is [20] 

m< 10^1^ eV. 

Standard cosmological solutions may be reproduced in 
massive gravity [18]. Further constraints on the param- 
eters of the model are imposed by the growth of pertur- 
bations and structure formation in the post-infiationary 
epoch [21]. 

In General Relativity (GR), a crucial role is played by 
the spherically symmetric vacuum solution to the Ein- 
stein equations — the Schwarzschild solution. This role 
is twofold: First, this solution describes the metric outside 
of spherical non-rotating bodies and gives rise, in the weak 
field limit, to the Newtonian gravity. It provides therefore 
a useful approximation in many astrophysical situations. 

Second, the Schwarzschild solution describes the result 
of a gravitational collapse, the black hole. Although the 
existence of black holes has not yet been directly con- 
firmed, there exists an indirect evidence that some of bi- 
nary stellar systems contain black holes as one of the com- 
panion [22, 23], and that many galaxies, including the 
Milky Way, harbor super-massive black holes in their cen- 
ters [24, 25]. It is conceivable that black holes will be 
directly observed in the near future, and that their prop- 
erties, including the metric configuration near the horizon, 
will be quantitatively tested [2G], thus providing a probe 
of GR in a fully non- linear regime. 



The Schwarzschild metric, together with properly ar- 
ranged scalar fields is a solution to Einstein equations in 
massive gravity as well [27]. However, the properties of 
black holes are, in general, expected to be different. In 
particular, rotating black holes are certainly modified, and, 
more generally, black holes are expected to have hair [27]. 
The possible existence of black hole hair in massive grav- 
ity models suggests that there might exist spherically sym- 
metric solutions other than the Schwarzschild one. In this 
paper we construct explicit examples of such solutions. 

We found a new class of vacuum spherically-symmetric 
solutions in massive gravity which depend, in addition to 
the mass M (equivalently, the Schwarzschild radius), on 
one more parameter which we call the "scalar charge" S. 
At zero value of this parameter the standard Schwarzschild 
solution is recovered, while at non-zero values of S the 
Schwarzschild metric gets modified. The modified solu- 
tion is non- linear at all distances; it cannot be obtained in 
the linear approximation. Similar phenomenon has been 
previously found in the context of bi-gravity models [28] . 
The new solutions may have event horizons and are, there- 
fore, candidates for modified black holes. We found both 
analytical and numerical examples of such modified black 
holes. 

The analytical solutions found in a particular massive 
gravity model show a variety of different behaviors. De- 
pending on the parameters of the model, the deviation of 
the metric from the Minkowski one may decay at infin- 
ity as 1/r or slower. The solutions thus may have finite 
or infinite ADM [29] mass, respectively. In the case of a 
finite mass, this mass may be positive or negative depend- 
ing on the sign of M . In either case the singularity at the 
origin may be hidden by the horizon. The solutions with 
negative mass exhibit an anti-gravitating behavior at large 
distances from the center. 

The modified gravitational potential may decay slower 
than 1/r in a certain distance range. Such a behavior 
would mimic the effect of the dark matter. Interestingly, 
the modification depends not only on the parameters of 
the model, but also on the scalar charge and the mass 
of the solution. Thus, the apparent amount of the "dark 
matter" may be different for objects of the same mass, in 
contrast with other models of modified gravity. 

This paper is organized as follows. In Sect. 2 we briefly 
review the massive gravity model and summarize previous 
results about the Newtonian potential in this model. In 
Sect. 3 we present the spherically symmetric ansatz and 
reduce the Einstein equations to four ordinary differential 
equations. In Sect. 4 we choose a particular model and 
find analytical solutions to the Einstein equations in this 
model. In Sect. 5 we demonstrate, by numerical compu- 
tation, that modified black hole solutions exist in more 
realistic massive gravity models which are attractors of 
the cosmological evolution. Finally, Sect. 6 contains the 
summary and discussion of our results. 



2. The massive gravity model 

In this paper we consider the massive gravity model de- 
scribed by the following action [IG], 

S = j Axy~g [~Ml{R + £,n + A^^] . (1) 

The first two terms are the curvature and the Lagrangian 
of the minimally coupled ordinary matter; they comprise 
the standard GR action. The third term describes four 
scalar fields (jP, (jf whose space-time dependent vacuum ex- 
pectation values break spontaneously the Lorentz symme- 
try. These fields are minimally coupled to gravity through 
a derivative coupling; they will be referred to as the Gold- 
stone fields. We consider the functions T which depend 
on two particular combinations of the derivatives of the 
Goldstone fields, T ^ T {X, W'^), where 

^ - A4 ' 

A4 A^X 

The constant A has the dimension of mass. The model is 
understood as the low-energy effective theory valid below 
the scale A. 

The vacuum configuration has the form 

g,.. = V^u, 0° = aA2t, <l>' = bA^x\ (2) 

where a and b are two constants determined by the re- 
quirement that the energy-momentum tensor associated 
with the four scalar fields vanishes in the Minkowski back- 
ground. The configuration (2) is, therefore, a solution to 
the Einstein equations. The constants a and b may be set 
to one by the redefinition of fields, which we assume to be 
the case in what follows. 

For functions J- which are invariant under rotations of 
the Goldstone fields cj)^ in the internal space (i.e., those de- 
pending on through three combinations Wn = Triy" , 
n = 1,2,3), the background (2) preserves the rotational 
symmetry. The Lorentz symmetry is, in general, broken. 

The action of the Goldstone fields with functions J- de- 
pending only on X and W^^ is invariant under the follow- 
ing symmetry, 

where ^re arbitrary functions of Because of this 
symmetry, the behavior of perturbations about the vac- 
uum (2) is non-pathological, i.e. there are neither ghost 
nor rapid instabilities [16]. The spectrum consist of two 
tensor modes (graviton polarizations) only, which are, in 
general, massive. The graviton mass scale is m ~ A^/Mpi. 

During the cosmological evolution, the Universe de- 
scribed by the action (1) is driven to an "attractor" point 
[18] 

T {X, W) T {Z'^) where Z'^ = X'^W'\ 
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7 being a constant. At the attractor point the theory 
possesses an additional symmetry, 



(3) 



Models described by the action (1) with the function = 
JF(Z'^) have been studied more intensively [18, 21, 30]. 
In particular, it has been shown that for —1 < 7 < and 
for 7=1 the cosmological perturbations in these models 
behave identically to those in GR [21]. For other values 
of 7 the behavior of the perturbations may or may not 
reproduce that of GR depending on the initial conditions. 

Another reason to study models characterized by the 
function F = T (Z*^) comes from the analysis of Newto- 
nian approximation. The gravitational potential of a static 
source in the model (1), in the linear approximation, has 
the form [18] 



where Gn = (STrMpj) is the Newton's constant and /i 
is a constant of order of the graviton mass whose value 
depends on the particular form of the function J- . This 
constant vanishes at the point where the symmetry (3) 
holds, that is, where T = T [Z''^). Thus, at the point of 
the attractor the Newtonian potential remains unmodified. 

It has been shown that the standard Schwarzschild met- 
ric is a solution to the Einstein equations in massive gravity 
models possessing the symmetry (3) [27], with the scalar 
fields given by 



,ln 



where Vg is the Schwarzschild radius of the black hole. 
On the other hand, the metric of a rotating black hole 
is necessarily modified [27] from its standard GR (Kerr) 
form. 



3. Static spherically symmetric ansatz and equa- 
tions 

A static spherically symmetric configuration in the mas- 
sive gravity model (1) can be written in the following form, 

ds2 = a{r)At^ + 25{r)AtAr - P{r)dr'^ - K(r)dr2^ 
<IP = K^[t + h{r)], 
A^x* 

<P' = 0(0—. 

r 

This field configuration is invariant under two residual co- 
ordinate transformations. The first one is an arbitrary 
change of the radial coordinate r ^ r' = r' (r) , which al- 
lows to set either k — r"^ or ip — r. The second one consist 
in redefining the time variable t t' = t + r (r). This last 
transformation allows to cancel either S (r) or h (r) . We 



choose the conditions k = r'^ and S = 0. Thus, we get the 
following ansatz. 



a{r)dt'^ - /3{r)dr^ 
A^ [t + h{r)], 

r 



sm' 



</>(r) 



(4) 



As compared to GR, this configuration contains two addi- 
tional radial functions h(r) and (f>{r). 

As has been pointed out in Sect. 2, the rotational in- 
variance of the vacuum (and likewise, of the ansatz (4)) 
requires that the function depends on W^^ through three 
combinations w„ = Tr(VF"). These combinations are ex- 
pressed in terms of the radial functions and their deriva- 
tives as follows, 

wi = - (/i + 2/2) , 

W3 =- {f! + 2/1) , 
where the two functions /i and /2 are 



and 



X = 



a(3X' 

(3 - ah'^ 
a(3 



/2 = 



In these expressions and in what follows, the prime denotes 
the derivative with respect to the radial coordinate r. 

After fixing the ansatz (4), the Einstein equations reduce 
to the following four equations: 



1 



= %\ 



where Q^i, is the Einstein tensor and is the energy- 
momentum tensor of the four Goldstone fields. The other 
six equations are identically satisfied. The explicit expres- 
sions for the components of Q^^ and 7^" are given in the 
Appendix. 

Consider first the equation T^*" = 0. Assuming h' 0, 
this equation gives 



(5) 



where = dT/dX and Ti = dJ-/dwi. Furthermore, the 
time and radial components of the energy-momentum ten- 
sor differ by the quantity proportional to eq. (5). There- 
fore, when this equation holds one has Tq = T^. This 
implies that Qq — or, equivalently, 

a(r)/3(r) = l. 
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in full analogy with the Schwarzschild solution in GR. 
Hence, the Einstein equations reduce to the following four 

equations, 



1 = 
= 

= 







a' 


a- 1 


1 

r 








a' 


a" 


1 

r 


~2 ^ 



m 
~2 



(T - 2XTx) 







~wiTi - 2W2T1 ~ SwaJ^s] , 

= XTx + /i (^1 - 2/1^2 + 3/2^3 
2 



(6) 
(7) 



(8) 

(9) 



where rr? — A^/Mpj. For a generic function JF, this system 
of equations is well defined. Indeed, since the function h{r) 
enters the equations (6-9) only through the variable X, 
one may consider X as an independent variable instead of 
h(r). Then the fourth equation allows to find (\) in terms 
of X, while the first equation gives /3 in terms of a. The 
second equation then gives X in terms of a and the third 
equation allows to determine a as a function of r. 

4. Analytical example 

Finding analytical solutions of the non-linear system of 
equations like (6-9) is impossible for a generic function T . 
So, in order to get some insight into the behavior of the 
solutions, let us choose the function T in such a way that 
the resulting equations are solvable analytically. 

Consider the function T of the following form, 

= co(^^ + w'i) (10) 

+Ci {w\ — 3u'iW2 — 6W1 + 2u'3 — 12) , 

where cq is an arbitrary dimensionless constant and c\ = 
±1 (the numerical value of c\ can be absorbed into the 
constant A). The coefficients inside the parentheses are 
chosen in such a way that the vacuum (2) is the solution 
to the Einstein equations at a = h ~ 1. Our example 
contains, therefore, a single continuous free parameter cq. 
Two additional constraints should be imposed on cq. The 
first one comes from the requirement that the graviton is 
non-tachyonic. This translates into the inequality 



Co - 6ci > 0. 



(11) 



The second condition is necessary to ensure that scalar 
modes with pathological behavior do not reappear upon 
addition of higher-derivative terms (recall that the model 
( 1 ) is understood as the low-energy effective theory, so such 
terms are generically present). This imposes certain con- 
straint on the function T and its derivatives (cf. eq. (70) 
of Ref. [16]), which in the case at hand implies 



Co > 0. 



(12) 



We should stress that the choice of the functional form 
of T is by no means unique. The particular form (10) has 
been chosen in order to simplify the solution of the field 
equations, as will become clear below. 



^.1. The static spherically symmetric solution 
Let us start with eq. (9) which takes the form 

1 



= 



X 



-Co + f • 



6ci 



Co - 6ci 



Because of our choice of the function J-, this equation con- 
tains the variable X as an overall factor only. Hence, it 
reduces to a closed differential equation for cj). The solution 
to this equation is 

(j) — br, 



where the constant b satisfies the equation 
0= (&2-l)(6fe* + 662 + co/ci). 



(13) 



We are interested in real positive values of b (the case 6 < 
can be reduced to 6 > by the inversion of coordinates) . 
If co/ci > there is only one such solution, 

b^l. 

If co/ci < there exists another one, 
1 / 



b = 



V2 



(-1 + v/l-2co/3ci) 



1/2 



Thus, at negative cq/ci we have two different branches of 
solutions. 

The remaining equations (7) and (8) can be written as 
follows. 






where 



a 
r 



1 



a" + \ 



3Ac -I- rn^co 1 



"^-3A. 
r 



1 

'x 

(2 + A), (14) 



A 



12^6 ^. 

Co' 

2„ fu& 



Ac = 2m^ci(6^-l). 

Recall that, according to our normalization, the constant 
Ci only takes two values, ci — ±1. 

Eq. (14) is a linear inhomogeneous equation for a. Its 
general solution can be found analytically. Making use of 
this solution and integrating the remaining equations one 
finally obtains: 



a (r) 
/3(r) 

h{r) 
<j,{r) 



1 



r r 
l/a(r). 



S 



± 

br. 



dr 
a 



1 



S A-1 



(15) 

n 1/2 



Here and S are two integration constants: is the 
Schwarzschild radius, while S* is a scalar charge whose 
presence refiects the modification of the gravitational in- 
teraction as compared to GR. At 5* = the solution (15) 
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reduces to the conventional Schwarzschild solution describ- 
ing a black hole of the mass M — r<,/(2G'Ar). 

The behavior of the metric at r ^ oo depends on the 
constants cq and ci = ±1. If ci = 1, we must also take 
Co > 6, and the only solution of eq. (13) in this case is 
6=1. Hence, we have Ac = 0, A < 0, and the metric is 
growing at infinity as Sr'^' . Such solutions do not describe 
asymptotically fiat space. 

If ci = —1, we must take cq > in order to satisfy 
eqs. (11) and (12). In this case A > 0. Two branches 
of solutions exist, one with zero, and one with non-zero 
cosmological constant, which can be positive or negative 
depending on the numerical value of cq. In what follows 
we consider the branch with Ac = (6 = 1). 

The behavior of the solutions is determined by the two 
integration constants M and S, and the value of the pa- 
rameter A. If A < 1, the third term on the r.h.s. of eq. (15) 
dominates at large distances. The ADM mass of these so- 
lutions is infinite. We do not discuss them further here. 

If A > 1, the standard Schwarzschild term dominates at 
infinity. The ADM mass of such solutions is equal to M. 
The solutions with positive (negative) M have attractive 
(repulsive) behavior at infinity. 

At the origin r = both terms proportional to M and 
S are singular, so the metric always possesses a singularity 
unless M = S = 0. This singularity may or may not be 
hidden by the horizon depending on the signs and values 
of M and S. The solutions possessing the horizon are 
candidates for modified black holes. 

The horizon is always present if both M and S are posi- 
tive. Such black holes have attractive gravitational poten- 
tial at all distances, which is stronger than for a conven- 
tional black hole of the mass M. The horizon size of the 
modified black hole is larger than — 2G]\[M. 

If M > and S <0, the presence of the horizon depends 
on the relative values of S and M. It exists for sufficiently 
small S. Defining, on dimensional grounds, the mass pa- 
rameter s associated with S by the relation \S\ = s~^, the 
existence of horizon requires that 



the horizon. The attraction changes to repulsion at 



sM > 



2G 



N 



A- 1 



A - 1 
A 



(16) 



The Newtonian potentials for solutions satisfying and not 
satisfying the condition (16) are shown in Fig. l-(a) and 
(b), respectively. When the horizon exists, the gravi- 
tational field is attractive all the way to the horizon. 
The attraction is weaker than in the case of the usual 
Schwarzschild black hole of mass M, and the horizon size 
is smaller. The behavior of the gravitational force with 
distance mimics that of the smaller-mass black hole plus 
a continuous distribution of "dark matter" , with the to- 
tal mass enclosed within the radius r approaching M as 
r ^ oo. 

Finally, at Af < and S > the modified black hole 
anti-gravitates at large distances and gravitates close to 



r = r^, — 



\S 



1 

A- 1 



The corresponding Newtonian potential is shown in Fig. 1- 
(c). 

A remark is in order at this point. In the conven- 
tional general relativity the constant or equivalently, 
the black hole mass M, is also a free parameter which 
can, in principle, be positive or negative. In GR, how- 
ever, only positive values make sense for the following rea- 
sons. First, negative-mass Schwarzschild solutions pos- 
sess naked singularity at the origin, which is physically 
unacceptable. Second, the conventional matter satisfies 
the null energy condition which ensures that any compact 
spherically-symmetric matter distribution has a positive 
mass [31]. None of these arguments go through in the case 
of massive gravity. Fig. l-(c) gives an example of solu- 
tion with repulsive behavior at large distances and with- 
out naked singularity: as for a conventional black hole, the 
singularity of this solution is hidden behind the horizon. 
The positivity of energy is also not expected in massive 
gravity. This is related to the fact that the background 
(2) breaks time translations, and only the combination of 
the time translations with the shifts of 0° by a constant re- 
mains unbroken. In this respect the massive gravity model 
is exactly analogous to the ghost condensate model [15], 
where the negative-energy states have been constructed 
explicitely [32]. 

4-. 2. Correspondence with linear analysis 

The solutions found in the previous section have the 
asymptotic behavior different from that obtained in the 
linear perturbation theory in Ref. [20]. In order to com- 
pare our results with those of Ref. [20] , let us discuss them 
in the gauge where h{r) — 0, 5{r) ^ (cf. Sect. 3). In this 
gauge the perturbation theory of Ref. [20] corresponds to 
assuming that the variations of all the metric components 
are of the same order. In other words, they are formally 
assigned a small parameter e to the first power. The solu- 
tions described in Ref. [20] satisfy the Einstein equations 
expanded to the linear order in e. 

The solution (15) is not of this type. Transforming it 
into the gauge h{r) = one finds that gtr — 5{r) does not 
decay as fast as the perturbations of other components, for 
instance, as a — 1. In fact, in the equations expanded in 
powers of perturbations the terms of order 5^ balance those 
linear in a — 1. In other words, in the formal expansion of 
the solution of Sect. 4 in powers of the small parameter e 
the perturbation 5 should be assigned the order rather 
than e. Hence, the solution (15) is non-linear even at large 
distances from the center. A similar phenomenon has been 
observed in the context of bi-gravity models in Ref. [28]. 

Another difference between the solution (15) and the 
solution to the linearized equations of Ref. [20] is that the 
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Figure 1: The deviation of goo from one (proportional to the Newtonian potential $ = (goo — l)/2) for tree different choices of the integration 
constants. Figs, (a) and (b) correspond to Af > and S < with the numerical values satisfying (a) and not satisfying (b) eq. (16)- Fig. (c) 
represents solution with M < and S > 0. 



former is static, while in the latter only metric components 
are static (in the gauge goi — 0)- The scalar fields have 
time dependence which may be viewed as an accretion of 
a fluid with zero energy-momentum tensor. 



4-.3. Gravitational field of a star 



In general relativity, one may relate the mass of a star 
to an integration constant of the vacuum solution in the 
exterior space by matching the interior and exterior solu- 
tions at the star surface (see, e.g., Ref. [33]). In massive 
gravity, one may try to use the same approach to deter- 
mine the scalar charge of an ordinary star. The analytical 
solution in the interior region is required for the matching 
procedure. 

The star is described, to a good approximation, by a 
diagonal energy-momentum tensor fj^ = (p, —p, —p, —p), 
where p and p are the energy density and pressure inside 
the star, respectively. This energy- momentum tensor is as- 
sumed to be responsible for the external gravitational field 
described by eq. (15). Since there is no direct coupling be- 
tween the ordinary matter and the Goldstone fields, fj^ 
must be conserved separately, V^iJ;^ — 0. For simplicity, 
we take the energy density to be constant at r < i?, where 
r = i? is the surface of the star, and zero outside. The pres- 
sure p cannot be chosen independently; it is determined by 
the conservation of tf,. 

Because of the spherical symmetry, the ansatz (4) holds. 
The Einstein equations in the interior of the star are ob- 
tained from eqs. (6-8) by adding the contributions of the 
energy-momentum of the star, while eq. (9) remains un- 
changed. The resulting set of equations can be solved an- 



alytically. The solution reads 

51 S2 , . 2 



a (r) 



1 



Arr' 



r 

p I r 



M 



pi 



R 



6 



')+o{p-)^ 



h{r) 
<j,{r) 



1 ^ + Acr - 



3M: 



± 



dr 

a 



pi 



1 I *2 A - 1 



n 1/2 



For simplicity, we have expanded the first equation in pow- 
ers of p, while the other relations are exact. Since the ge- 
ometry inside the star is regular, the integration constants 
si and S2 must be set to zero. 

The interior solution has to be matched with the solu- 
tion (15) at r = i?. It is convenient to match the variable 
X which equals 1 in the interior region. In the gauge 
h{r) — this variable is nothing but the g^^ component 
of the metric. Hence, it must be continuous. Making use 
of eqs. (15) one can see that the continuity oi X at r = R 
requires that S = Q. Therefore, the scalar charge of an 
ordinary star is zero. 

It remains an open question how objects (e.g., black 
holes) with 5* ^ can be created. The argument given 
above does not apply to time-dependent configurations, so 
it is possible that a non-zero scalar charge may be acquired 
during the gravitational collapse. 



5. :F(Z^J) models 

As mentioned earlier, models characterized by the func- 
tion of a single variable Z^^ — X^W^^ are of a particular 
interest. We discuss in this section the exact static spher- 
ically symmetric solutions in these models. Our goal is to 
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demonstrate that the solutions found earher are not spe- 
cific to the particular form of the action (10) and exist also 
in the models obeying the symmetry (3). 

In Sect. 4 the analytical solutions of eqs. (6-9) were 
obtained by choosing the function T in such a way that 
the dependence on X factors out in eq. (9). Since now T 
has only one argument, the derivatives of T with respect 
to X and W^^ are no more independent. For this reason we 
did not succeed in constructing non-trivial examples where 
the Einstein equations are solvable analytically. Hence, 
to demonstrate the existence of solutions we have to use 
numerical methods. 

Consider the following function T , 

z\ - 6ziZ2 + 82:3' 
3 



goo - 1 



T = Co hzi + 2 + 



(17) 



2ci zf - 2z2 - 1 - 2- 



6ziZ2 + 8Z3 



were z„ = Tr (Z" ) are tree independent scalars made out of 
Z*^ , Zn = {X'^Y^ Wn- The coefficients in front of individual 
terms have been adjusted so that the flat metric and the 
scalar fields given by eq. (2) solve the field equations at 
a = b — 1. For the vacuum configuration (2) one has 
Z*-' — —S"^^ . We are interested in solutions to the field 
equations that asymptote to this vacuum state. 

In addition to the adjustments already made, the fol- 
lowing inequality should be imposed on the coefficients Cq 
and ci to ensure that the graviton is non-tachyonic. 



Co - 2ci > 0. 



(18) 



This guarantees that the square of the graviton mass is 
non-negative. Moreover, this inequality is sufficient for 
the absence of pathological scalar modes which may ap- 
pear upon addition of higher-derivative terms. As in the 
previous example, the overall scale of the coefficients cq 
and ci can be absorbed in the parameter A, so without 
loss of generality we may set ci = ±1. 

For this class of models, the field equations (G-9) may 
be viewed as equations for a{r), f3{r), ^(r) = X"' fi — 1 
and ip{r) = X^ f2 — 1. Then eq. (6) gives /? in terms of 
a, while eq. (9) enables to express ip in terms of ^. The 
two remaining equations form a coupled set of non-linear 
equations for a and 1^ ; they have to be solved numerically. 

The numerical solutions are shown in Fig. 2 for different 
value of the parameters cq and ci. For all these graphs, 
we have assumed that the external horizon is located at 
r = 1 and that ^ — 100 at the horizon. The large value 
of ^ is chosen in order to make the difference between the 
modified solution and the Schwarzschild solution visible 
on the plot (large values of ^ correspond to large scalar 
charge S of the previous section). 

The plots show the behavior qualitatively similar to that 
discussed in Section 4. In particular: 

• if Co = Ci = — 1, the Newtonian potential 2<i> = a — 1 
is attractive at short distances and becomes repulsive 
at larger distances; 




- Schwarzschild Solution 
-c„ = - 1,C| = - 1,7 = -l,m = 0.001 
_ e„ = 3, t, = 1 , 7 = -1, m = 0.001 



Figure 2: The deviation of goo from one for three different cases; 
the usual Schwarzschild solution (solid line) and two solutions cor- 
responding to different values of the parameters of the function (17) 
(long-dashed and short-dashed lines). The integration constants of 
these solutions have been chosen such that the external horizon is 
located at r = 1. 



• if Co = 3 and ci = 1, the Newtonian potential is at- 
tractive outside of the horizon, and becomes repulsive 
close to the singularity. 

The deviations from the Schwarzschild metric are larger 
for larger values of the "scalar charge" (parameterized by 
the value of ^ at the horizon) . The Schwarzschild solution 
is recovered at ^ 0. 

6. Discussion 

To summarize, there exist spherically symmetric vac- 
uum solutions in massive gravity models which depend on 
two integration constants, the mass M and an extra pa- 
rameter S which can be called the "scalar charge". At 
zero value of the scalar charge the standard Schwarzschild 
solution is recovered, while at non-zero S the metric is 
modified with respect to the Schwarzschild case. 

The solutions having non-zero scalar charge exhibit 
much reachcr behavior than the Schwarzschild solution in 
GR. As can be seen from the explicit example of Sect. 4, 
both the short and long distance behavior may be modified 
depending on the parameters of the model. 

Unlike in General Relativity, the solutions may have a 
negative ADM mass. Such solutions have repulsive gravi- 
tational interaction at large distances. At short distances 
the repulsion may change to attraction and give rise to 
the horizon, hiding the singularity at the origin. Such so- 
lutions represent anti-gravitating black holes. 

In the case of a positive ADM mass, the S'-dependent 
contributions may make the gravitational attraction 
weaker at short distances (cf. fig. l-(a)). In this case 
the gravitational force decays with distance slower than 
l/r^, thus mimicking the presence of dark matter. Inter- 
estingly, solutions with the same value of M but different 
scalar charge S have different behavior, which corresponds 
to different amount of the apparent "dark matter" . This 
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is in contrast with other models possessing modifications 
of the gravitational potential [4, 6, 10], where the mod- 
ification of the gravitational force is determined by the 
parameters of the model. 

It is currently an open question how objects with non- 
zero scalar charge may be created. As has been argued 
in Sect. 4.3, the absence of direct couphng between the 
Goldstone fields and ordinary matter results in zero scalar 
charges of static matter distributions. Thus, the gravita- 
tional field of ordinary stars is described by the S = 
solutions, i.e., by the standard Schwarzschild metric. This 
may be not the case for black holes, especially the super- 
massive black holes in the centers of galaxies, which may 
be of primordial origin [34, 35]. In any case, this question 
requires further investigation. 

Another open question is the stability of the modified 
black hole solutions. Several kinds of instabilities may be 
present. Among perturbations of the solutions there may 
exist unstable modes with the characteristic time scale of 
order of the horizon size; in this case the interpretation 
in terms of black holes is not possible. Second poten- 
tial source of problems is generic presence of the higher- 
derivative terms not included in the action (1). One has to 
check that the 2-parameter family of modified black holes 
survives their inclusion. By analogy with the ghost con- 
densate case, one may expect that these terms produce at 
least a slow Jeans-type instability [15], which is not, how- 
ever, dangerous for the black hole interpretation. Finally, 
the presence of negative mass solutions may lead to insta- 
bilities of the quantum- mechanical nature similar to those 
found in Ref.[36]. 

To conclude this list, let us mention also the solutions 
satisfying h = which were not considered in this paper. 
In this case the ten Einstein equations reduce to three 
equations for a, (3 and 4> which form a (generically) well- 
defined system. It remains to be seen whether this sys- 
tem has asymptotically flat solutions. In any case, the 
Schwarzschild solution does not belong to this class which 
is characterized by aP ^ 1. 
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Appendix 

For the ansatz (4) discussed in Section 3, the non-zero 
components of the Einstein tensor are given by the follow- 



ing relations 



r 

a + ra' 
a(3 



1 

4r 



a' + ra" 
al3 



2a + ra'\ ' 



The energy-momentum tensor of the Goldstone fields is 
expressed through 



1 



X 



where Tx = dT jdX and 



A4 



For functions T which are invariant under rotations of the 
Goldstone fields internal space, the derivatives of T with 
respect to W^^ are given by 

Fit 

where Ti = dT jdwi. Therefore, the components of the 
energy-momentum tensor which are not identically zero 
are given by 



A^ 



A^ 



1 



2 

1 dT 



T-[Tx 



dT V'V^ 
dW^i A2 

dT V'V^\ h'^ 



p dW^i \ X 



7-v ^ _iA4jr- 
f 2 



dT 



d'^ch-'d. 
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